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The p and A mesons in a strong abelian magnetic field in SU(2) lattice gauge theory. 

E.V. Luschevskaya^ and O.V. Larina^ 

^ ITEP, B. Cheremushkmskaya 25, Moscow, 117218 Russia 
(Dated: March 26, 2012) 

We calculate the correlators of vector, axial and pseudoscalar currents in the background of 
a strong abelian magnetic field in a pure SU{2) gauge theory with improved action. Using the 
maximal entropy method and fitting the correlators by a coshinus function we extract masses of p 
and a meson ground states in different spatial directions. We found that with the growth of the 
external magnetic field the p meson mass along the magnetic field P33 decreases, while the masses 
of the components perpendicular to the magnetic field, P22 and pn slowly increase. The axial 
meson mass along the magnetic field also diminishes with the value of the field. However in the 
perpendicular directions to the magnetic field the mass components an and 022 do not change up 
to magnetic field strength ~ 1 GeV and stay equal each other within the error bars. 



PACS numbers: 11.30.Rd; 12.38.Gc; 13.40.-f 



I. INTRODUCTION 

Quantum Chromodymanics in the background of an 
external magnetic field is a rich realm for explorations 
as it presents a huge interest for physicists because of 
many reasons. At first the strong magnetic fields can 
be related to the formation of the early Universe. It is 
not clear whether these fields were fundamental or occur 
as a result of a mutual motion of charged particles at 
that era. It is suggested, that the magnetic fields ^ 2 
GeV existed in the Early Universe during the electroweak 
phase transition [1]. It is known that there are objects 
in the cosmic space, called magnetars or neutron stars, 
having the magnetic field in their cores of the order ^ 1 
MeV. 

Many amazing effects were discovered experimentally 
and approved theoretically in the recent time. Experi- 
ments by the STAR collaboration at RHIC in An- An col- 
lisions showed that the chiral magnetic effect takes place 
in the noncentral heavy- ion collisions 043 ■ This effect 
was confirmed in ALICE experiment at LHC. The strenth 
of the magnetic field can reach Ibrn^ ^ 290 MeV"^ [a|, 
that is of the order of hadronic scale. The magnetic field 
appears because of the mutual motion of the ions (it is 
known that a moving charged particle creates an electric 
and magnetic field around itself) and creates an asymme- 
try of charged particles emitted. The essence of this ef- 
fect is an appearance of a charged asymmetry of charged 
particles, emitted from the different sides of the reaction 
plane. So, it occurs that such a strong magnetic field can 
be created in the Earth laboratories and it is possible 
to explore the quark - hadronic matter in such extreme 
conditions, thus making this system very interesting. 

A strong enough magnetic field can lead to the mod- 
ification of the QCD phase diagram. Phenomenological 
calculations point out that the temperature Tc of the 
transition from the confinement phase to the deconfine- 
ment one increases with the value of B and the transition 
becomes of the first order Q . An increase of Tc is also 
predicted by models of Nambu-Jona-Lasinio type: NJL, 



EPNJL, PNJL and PNJLg i, Gross-Neveu model 
[3, [Ifll ■ The first prediction from lattice simulations for 
two flavour QCD about the behavior of the deconfining 
temperature and the chiral symmetry restoration with 
increasing of the magnetic field was made in |ll| . Con- 
trary to this result the lattice simulations with Nf = 2 + 1 
fiavours [l^l shows the decreasing of Tc with magnetic 
field. However, the chiral perturbation theory predicts 
the opposite effect of decreas ing the transition tempera- 
ture with the magnetic field jl3{. 

The same have also been studied a few years earlier, 
but for the case of a chromomagnetic (instead of mag- 
netic) field, in [ll{l3. 

The external magnetic fields lead to the enhancement 
of the chiral symmentry breaking in the background of 
the field [l342l|. The analytical calculations (2342i| 
predict a linear increase of chiral condensate with the 
magnetic field in the leading order. Lattice simulations 
show that the chiral condensate depends on the strength 
of the applied field as a power function with exponent n = 
1.6 ± 0.2 [Ml. AdS/CFT approach shows the quadratic 
behaviour [26|. 

In |27| was provided the first numerical study of mag- 
netic catalysis in QCD with dynamical fermions. 

Recently it was obtained from the Nambu-Jona- 
Lasinio model that in strong magnetic field the vacuum 
becomes a superconductor '28i| . The transition to the su- 
perconducting state occurs at the critical value of mag- 
netic field equal to Be = m'^/e ~ 10^^ Tesla. The super- 
conductivity is realized along the axis of the magnetic 
field only, and is accompanied by the condensation of 
p mesons. Further investigations in this direction have 
been done in lattice investigations and [2^, |30| and point 
out at the presence of a superconducting phase at strong 
magnetic fields. In our present work we explore the be- 
haviour of the p and a meson ground states in the abelian 
strong magnetic field. We observed a deacrese of the me- 
son masses in the direction of the magnetic field. We 
work in SU(2) gluodynamics because we are interested 
in qualitative predictions for the beginning. This model 
contains the main properties of the full QCD theory and 



easy for numerical simulations. 

The paper is organized in the following way. In Sec- 
tion |ll] we describe the technical and numerical setup of 
our simulations. We discuss measured observables in Sec- 
tion mil Section IIVI is devoted to the methods which we 
use. The results of our calculations are presented in Sec- 
tion lVIl Our results are compared with the previous ones. 



II. TECHNICAL DETAILS AND NUMERICAL 
SETUP 

For the generation of SU{2) quenched configurations 
we use the tadpole improved Symanzik action 
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where Spi^rt ~ (l/2)Tr(l — Upi^rt) is the plaquette (de- 
noted by pi) or 1x2 rectangular loop term (rt), uo = 
(W'lxi)^/^ = ((l/2)TrC/p,)i/4 is the input tadpole fac- 
tor, computed at zero temperature [31|. This action sup- 
presses the ultraviolet dislocations, leading to unphysical 
near-zero modes of the Wilson-Dirac operator. 

For the calculation of the fermion spectrum we use the 
chiral invariant massless Neuberger operator [32[. It is 
a nontrivial construction, allowing to preserve the chiral 
symmetry on the lattice and exploit the perfect chiral 
properties of the Dirac eigenmodes. 
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where Dw = M—p/a is the Wilson-Dirac operator with a 
negative mass term p/a, a is the lattice spacing, M is the 
Wilson hopping term with r = 1. Anti-periodic boundary 
conditions in time and periodic ones in spatial directions 
are imposed to the fermionic fields in the overlap Dirac 
operator. Using the minmax polynomial approximation, 
we calculate the sign function. 



Dw/^/DWw = 75sign {Hw), 



(3) 



where Hw = 75 -Dm/ is the hermitian Wilson-Dirac op- 
erator. We treated the 30-50 lowest Wilson-Dirac eigen- 
modes separately for the calculation of the sign function. 
In the continuum the analog of this operator is D — 
I'^idfi — ^^/i)i satisfying the Dirac equation 



Dil^k = iAfeVfe- 
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Here A^ is the gauge field, Vfc is the eigenmode and Afe is 
the corresponding eigenvalue. The overlap operator Dov 
measures the spectrum of the eigenmodes and eigenval- 
ues of the probe quark, which we put into the gauge 
background {A^) of gluodynamics. The eigenmodes of 
the Dirac operator allow to construct operators and cor- 
relators, which we talk about in Section 3. 

We include the external uniform magnetic abelian field 
in the simulations. The abelian field interacts only with 



fermionic fields, so we introduce it only into Dirac oper- 
ator without any modification of the gauge action by the 
following exchange 
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Expression ([5]) is valid in an infinite space, therefore to 
combine it with the periodic boundary conditions on the 
lattice we have to introduce an additional z-dependent 
twist for fermions, as expalained in [33| . 

The magnetic field in our calculations is directed along 
the third axes, it's value is discretized. 
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because the lattice represents a finite 4 dimentional box, 
we have the only one type of quarks and take its charge 
q = — l/3e, which corresponds to the d-quark. The 
quantization condition posses a restriction to the mini- 
mal value of the magnetic field considering on the lattice. 
For example on the 16* with spacing a — 0.0998 fm the 
minimal value yeB — 535.64 MeV For inversion, we use 
a Gaussian source with radius r = 1.0 in lattice units in 
both spatial and time directions and a point sink (that 
is, quark position is smeared over a Gaussian profile). 

We performed the calculations at 14'* and 16^ lat- 
tices for different lattice cutoff, a = 0.0681, 0.0998 and 
0.138348 fm in order to control the effects of finite lattice 
spacing and volume. We used 300 configurations for the 
14** lattice with a = 0.138348 fm. However for the other 
lattices the statistics is small ~ 100 — 150 configurations 
and we are increasing it now to make the extrapolations 
later. We also perform the simulations at different quark 
bare masses from rriga = 0.01 till 0.5 in lattice units. For 
quark masses larger than 0.1 — 0.2 the procedure works 
well and we get the good results with the reasonable er- 
rors. 



III. OBSERVABLES 

We calculate on the lattice the observables of the type 

{^^{x)OMx)^Hy)02^{y))A, (8) 

where Oi , O2 = 75, 75 7/^, 7/^. In the Euclidean space 
"0^ = "0 |34| . These current-current correlators in the me- 
son channel can be represented in terms of Dirac propa- 
gators in the background of the gluonic fields and abelian 
constant magnetic field. In the continuum field theory (|8]) 
can be written in the following way 
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The first term in the difi'erence is the disconnected part, 
while the second one is connected . The disconnected 
part is very intricate for the calculations and has big 
errors. We check that the disconnected contribution is 
significantly smaller than the connected one and can be 
excepted from the consideration. So, we calculate only 
the connected part of the correlator ([9]) . 

We calculate the Dirac propagator on the lattice by 
inverting the massive Dirac operator 1/{D + m) in the 
subspace of M (30, 50 in our case) eigenvectors with the 
Xk lowest eigenvalues. It has the form 
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Taking into account pUj) one can reformulate © on the 
lattice 
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The first term in[5]is not essential as we have said above. 
The magnetic and gluonic fields are incorporated into the 
consideration in the way described above. 

The mass of neutral p-meson was extracted from the 
vector correlator of the vector currents {j^ {x)j^ {y))A, 
where jY {x) = 4'''{x)jfj.'4'{x), the mass of the neu- 
tral axial meson - from the correlator of axial currents 
{jt^{x)j^{y))A, where j^{x) = ■ip\x)-ff,j5'ijj{x) is the ax- 
ial current. The correlator {j^ {x)j^ (y)) a corresponds 
to TT-meson, here j^^ = ip'' {x)j^'ip{x) is pseudoscalar cur- 
rent. 



IV. METHODS 

We use two methods for the calculations of the masses. 
In the first method we consider the imaginary parts of the 
corresponding correlators and project the sink operator 
to zero spatial momentum. We get 



(V'HO, nOOi^(0, ^0^-^(0, 0)02^(0, 0)) A 
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For a finite lattice the values of Ek are discrete. For 
the single particle intermediate states the mass ruk — E^ , 
i.e. coincides with the low-lying energy value. We are 
interested in the ground state energyes of the mesons 
and its dependence from the magnetic field. 
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The first term 
the main contribution to the sum for the large values 
of lattice time slice nt. For this reason we fit the corre- 
sponding imaginary part of the correlators in the min- 
imum by the coshinus function of the form C{nt) = 
Acosh{meff{Tit — NT/2)) on the interval 4 < n^ < A^t~4 
for different nt and making averaging (fitting with differ- 
ent nt reduces the errors), Extract the mass ni^ff from 
the fit, Nt is the lattice temporal extention and get suf- 
ficiently good fits for the correlators within the range of 
error bars. 

The second method which we use is the Maximal En- 
tropy Method (MEM) ^. The imaginary-time Eu- 
clidean correlator G{t,p) = J (f x{O{T,x)O^0,6))e-'P^ 
is related to the spectral function p{uj,'p) according to 
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G{r,p)=J ^K{T,co)p{io,p). 



(13) 



In general p{uj,p) carries the all properties of mesons 
and hadrons. We consider p — and drope the p depen- 
dence in the following. The peak in the spectral function 
corresponds to the energy of the ground state. The kernel 
in P^ is given by 



K{T,i0) 
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where T is the temperature, r is the euclidean time, tu is 
the frequency. To extract the spectral function from p^ 
one can try to make an inversion. 

On the lattice this problem is ill-defined, because we 
calculate the correlator G{t) numerically at the discrete 
set of points r^ = Tmm + (i — l)a, where a is the lattice 
spacing, i — l,...Nr, Nr = l/{aT) is typically of the 
order of 0(10). Then we approximate the integral by a 
discrete sum at the points a;„ = nAa;, n — 1,..., N^^ and 
7V<^ is usually ~ 0(10^). We cutoff the integral ^ at 
some ijJmax- All the same the inversion becomes impos- 
sible. But the ideas of Bayesian probability theory allow 
to overcome this difficulty. 

The most probable spectral function p{u)) can be con- 
structed if we find the maximum of the conditional proba- 
bility P[p\DHam], where D is the data,7J is our hypoth- 
esis, a is a real and positive parameter, to, — m{uj) is the 
default model which we use. This procedure is equivalent 
to a maximization of the free energy F = L — aS, where 
S is the Shannon entropy term, defined by the following 
way 
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L is the standard likelihood function, the detailed ex- 
planation how to make a corresponding discretization on 
the lattice is given in [3a| • The parameter a balances the 
relative importance of the data and the prior hypothesis. 
We take a S [amim otmax] and average the data over this 
interval. This interval was choosen in such way that the 
results vary slightly, for approximately 10%. 



The kernel p^ contains divergence at a; = leading 
to unstable behaviour of the procedure at small energies. 
The Bryan's key idea was to redefine the kernel and spec- 
tral function 

K{uJ, t) - ^KiLo, t), piLo) = —p{uj) (16) 

so that K{uj,t)p{uj) = K{oj,t)p{uj) and apply the SVD 
theorem to the modified discretized kernel K{uJn,Ti), 
see [36|. We use this modified algorithm to determine 
the spectral function in the form 
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The column vectors Ui, (i = 1, .., A'^) are normalized 



N^ 



{Ui\Uj) = ^ Ui{Un)Uj{Un) = S^. 



(18) 



Ci are the coefficients and we set K{0,t) — 1. 

Therefore to reconstruct the spectral function p(w) we 
have to choose the default model 771(0;) in a correct way. 
The default model should describe correctly the high and 
low energy behaviour of the spectral function. Following 
the analysis of [33] we choose it in the form 



m{u)) — vriaUJ + rrih, vria 



GiNr/2) 



y2 



mb^ an: 



_3_ 

(19) 
1 for scalar and pseudoscalar channel, qh = 2 for 
vector and axial vector channel [38] . We also try another 
default model (constant function, ^ w^, vary the iria and 
nib), but the choice ([T9|) gives the best convergence of 
MEM. 



an 



C5 

CM 

E 




100 



FIG. 1: The squared mass of the neutral vr meson extracted 
from the C^ (nt) versus the bare lattice quark mass for 
14* lattice volume, lattice spacing 0.1383 fm, /3 = 3.1000 and 
different values of an external magnetic field. 
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FIG. 2: The same that in Figure[T] but for the lattice spacing 
0.1155 fm, /3 = 3.2000. 



V. VECTOR MESON MASS AT B = 

At this section we present the calculation of p 
meson mass in SU(2) gluodynamics for zero mag- 
netic field. We compare our results with the previ- 
ous one to be sure that our method works correctly. 
We measure the correlator of pseudoscalar currents 
C/^P^int) = (j^^(0,nt)j^^(0,0))^, where j^^(0,nt) = 
■0(0, nt)75^(0, nt) and calculate the mass of the lowest 
energy state of the neutral n meson for different quark 
masses, volumes and lattice spacings, see Figure [1] and [2] 
We get the linear dependence of the squared vr meson 
mass from the bare quark mass. It is known that the Chi- 
ral Perturbation Theory predicts such linear dependence 
according to the relation 



/>2 ^ ^-"(^V-), 



(20) 



where /^ is the pion decay constant, (ipip) is the chiral 
condensate, mj'^" is the renormalized quark mass. In 



the limit of zero quark mass we get massless pions. On 
the lattice we work with a bare quark mass. We see 
from the Figures [1] and [5] that the corresponding lines are 
shifted with respect to the origin of the coordinates and 
the value of quark mass renormalization slightly varies for 
the different lattice volumes and spacings within — 10 
MeV and its precise determination is impossible for small 
volumes. 

The mass of neutral p meson is calcu- 
lated from the correlator of vector currents 
C'^int) ^Jjj:{6,nt)j^{6,0))A, where j^{6,nt) = 
"0(0, ri()7^i/;(0, nt) from the coshinus fits Then we per- 
form the extrapolation to the infinite physical volume. 
We ignore the difference in the value of quark mass 
renormalization for different lattice volumes for a given 
lattice spacing and just take the averaged value. Later 
we see that the p meson mass is not very sensitive 
to the quark mass renormalization, at least the shift 
in quark mass for 10 — 20 MeV doesn't change nip. 
The extrapolation of p meson mass to the infinite 



physical volume limit was performed for several values 
of quark mass and two lattice spacings 0.1383 fm and 
0.1155 fm (Figures |4] and E]). The best fit gives the 
^/yphys volume dependence, Vphys = (clL)^, the power 
functions of the physical volumes give the worse fits. 
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FIG. 5: The mass of p meson for the different quark masses 
and two lattice spacings. The extrapolated values of masses 
are depicted by stars and corresponds to the physical quark 
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FIG. 3: The mass of the neutral p meson extracted from the 
C (nt) for several bare quark masses, lattice volumes and 
lattice spacing 0.f383 fm, /3 = 3.1000. Magnetic field equals 
zero. Stars correspond to the extrapolated values of mass. 



VI. THE MASSES OF MESONS AT B / 
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FIG. 4: The same that in Figure|3] but for the lattice spacing 
O.f f55 fm, /3 = 3.2000. 



The second step is the quark mass extrapolation. The 
corresponding dependences are shown at Figure [5] We 
take the value of rrip meson mass at the physical value 
of quark mass 3 MeV. The corresponding pseudoscalar 
mass ttItt — 150 — 250 MeV (in depence on lattice 
parameters). The vector meson mass rrip ~ 932 ± 12 
MeV for a = 0.1338 fm and 937 ± 12 MeV a = 0.1155 
fm lattice, that is in agreement with the previous 
evaluations. The corresponding points are depicted by 
stars. The extrapolations were performed to the quark 
mass mP^ = 3 MeV, taking into account quark mass 
renormalization on the lattice. In our consideration the 
light quarks are degenerate. 
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FIG. 6: The squared mass of the neutral tt meson extracted 
from the C^ {nt) versus the bare lattice quark mass for 
14* lattice volume, lattice spacing 0.1383 fm, j3 = 3.1000 and 
different values of an external magnetic field. 



At Figure |6] the squared pion mass is depicted for dif- 
ferent bare quark mass and different values of magnetic 
fields H = \feB. The results obtained from the Maximal 
Entropy Method agree with the ones from the coshinus 
fitting. Both methods give the linear dependence, that 
is legitimate for the the Chiral Perturbation Theory, be- 
cause this behaviour is in accordance with the relation 

From the (pn|) dependence and Figure [S] we conclude 
that the quark mass on the lattice in the background of 
the gluon and abelian magnetic field have to be renormal- 
ized and this renormalization depends on the magnetic 
field value (because the line is shifted for different values 
of B with respect to the origin of coordinates). 

We show the dependence of the tt - meson mass from 
the value of the squared magnetic field, see Figure [T] It 
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FIG. 7: The mass of the neutral tt meson extracted from 
the C^ (nt) versus the squared value of magnetic field for 
renormalized and nonrenormalized quark mass. 
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FIG. 9: The P33 mass along the direction of B, extracted from 
the C33^(nt) correlator for 14* and 16* lattices and different 
lattice spacings, obtained from Maximal Entropy Method. 
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FIG. 8: The P33 meson mass component along the direction 
of B, extracted from the C^^ (nt), versus the squared value 
of the field for 14* lattice volume, lattice spacing 0.1383 fm, 
P — 3.1000 for renormalized quark mass and for the bare one 
niqa = 0.01. The data are calculated by the MEM procedure 
and cosh-function fitting. 



appears that the mass of n is very sensitive to the quark 
mass renormalization, doesn't depend on the lattice vol- 
ume. For the magnetic fields less than 1 GeV we get 
the linear mass dependence from the magnetic field, the 
slope is negative, that is in accordance with the results 
of A.Smilga, obtained in chiral perturbation theory |24| . 
The value of the slope is smaller and differs from ChPT, 
because we explore the simpler SU{2) gauge theory with- 
out dinamical quarks and only qualitative behaviour has 
sence. 

The errors of MEM are very difficult to estimate, for 
the simplicity we take an ansamble of 0(10) results of 
MEM procedure and calculate an average and standard 
deviation for the best set of MEM parameters. Errors are 
calculated taking into account also the w-discretization. 

We calculate the components of the correlator of vec- 



FIG. 10: The same that at Figure |9] but from the cosh- 
function fit of the correlators. 



tor currents, the diagonal components are essentially 
nonzero, while the nondiagonal ones are zero within the 
error bars. The external magnetic field is directed along 
the third coordinate axes. The spatial diagonal com- 
ponents of the correlators perpendicular to the magnetic 
field are CY/{nt) = {jY {0,nt)jY (0,0)} a and 0^2""^ = 
UY{0,nt)jY{0,0))A, where jY{0,nt) = ^^(0, 0)7iV(0,nt) 
and so on. The component along the direction of the 
magnetic field is CY/ (rit) = {jY {0,nt)jY {0:0))a- We fit 
the correlators by the coshinus function and also apply 
the MEM for the extracting the spectral functions for 
the each correlator. The position of a peak at w 7^ of 
a spectral function corresponds to the mass of the state 
in our case. The Euclidean time CYq^ component is con- 
stant, in accordance with the conservation law of an elec- 
tric charge, d^jY = 0. The Cqo^ componant gives just a 
peak at a; = 0, we do not consider it here. 

Let's denote the masses extracted from the correlators 
perpendicular to the magnetic field by pn and P22 and 
the mass component along the field by P33. Figure [5] 
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FIG. 11: The pn meson mass extracted from the CYi^ (rit) 
correlator perpendicular to the direction of B for 14* and 
W^ lattices and for different lattice spacings from coshinus 
function fit the correlators. 



FIG. 13: The 033 mass, extracted from the C33 (nt) correlator 
along the direction of eB for 14* and 16* lattices and for 
different lattice spacings by means of MEM procedure. 
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FIG. 14: The 033 mass, extracted from the C33 {nt) correla- 
tor along the direction of B for 14* and 16* lattices and for 
different lattice spacings from MEM procedure. 



shows that the mass of P33, obtained from the C^^ {nt), 
along the direction of the magnetic field is insensitive to 
the renorm.alization. 

We see that the P33 masses obtained with Maximum. 
Entropy Procedure is a little bit lower that the masses 
extracted from the coshinus fit, it was to be expected 
because the MEM adjusts the correlators better than the 
coshinus fit does. Because the p mass doesn't depend on 
the renormalization we explore its behaviour for the bare 
quark mass mqa — 0.01. The fit at Figure |8] is a linear 
function in eB for the fields till 1 GeV . It's slope is 
equal to c ~ —0.5. In NJL approach [3^ the neutral 
p meson doesn't depend on the magnetic field and the 
slope ci = 1, for the charged p± the slope C2 ^ —0.6 for 
the interval [0, 0.4] GeV. 

At Figure [S] we depict the P33 mass, calculated from the 
{nt) correlator, but in another scale, versus veS, 
and for the larger values of the magnetic field. The same 



r'VV 
^33 



dependence, but obtained with coshinus fit the correla- 
tors, is shown at Figure [TOl We see that the p^z mass 
is decreasing with the growth of magnetic field and this 
takes place for the all lattices which we consider. How- 
ever near the 1 GeV the squared behaviour in B changes 
to the asymptotic one. We need the bigger lattices and 
statistics to perform extrapolations and check this be- 
haviour. 

The pii and P22 masses, obtained from the Cn and 
C22 correlators are shown at Figure [11] and [121 In con- 
trast to P33, obtained from C33, pn and P22 mass com- 
ponents increase with the field slowly. We observe the 
mass splitting because the external magnetic field create 
the anysotropy. In the absense of the magnetic field the 
masses pn, P22 and P33, extracted from the all correla- 
tors, coincide within the statistical errors. 

We also calculate the axial correlator for the different 
lattice parameters and magnetic fields and extract the 
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FIG. 15: The an, i — 1,2 masses, extracted from the Ca (ut) 
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lowest energy state. The 033 extracted from the C^2^{nt) 
along the direction of the magnetic field is depicted at 
Figure [T3](MEM) and [Til fcoshinus fit the correlators). 
We observe the decreasing of a33 mass component along 
the magnetic field. The peak around 1 GeV can be just 
a lattice artefact. 

In the perpendicular direction to the magnetic field 
an and 022 masses do not reveal any dependence on the 
magnetic field till 1 GeV. We can see this from Figure [T5] 
an, i = 1,2,3 are also independent on the quark mass 
renormalization. 



VII. LATTICE EXTRAPOLATIONS AT B 7^ 

On the lattice we work with discrete values of the exter- 
nal magnetic field. The value of magnetic field is deter- 
mined according to [7] and as it follows from this formula 
depends on the lattice volume (aL)^ and lattice spacing 
a. For this reason it's difficult to make lattice extrapola- 
tions to the continuum. At Figure [!§] the dependences of 
the ?Tip33 component (along the magnetic field) from the 
bare quark lattice mass are depicted to demostrate how 
we make a mass extrapolation. The third p mass compo- 
nent was shown here for the lattice volume 14'', the mag- 
netic field H = 6 in lattice units, but for different lattice 
spacings 0.0998, 0.1155 and 0.1383 fm. Cause the lattice 
spacings are different the values of magnetic field also dif- 
fer. We fit the dependences by mp = ao + ainiq + 02™^ 
function, which gives the best fit. We make a renormal- 
ization of lattice quark mass and after that extrapolate 
the mass of p meson to the point, where the physical 
quark mass is equal 3 MeV (marked by stars). As we 
mentioned above this value gives the mass of n meson 
close to the physical one. 

The p meson mass for the larger lattice volume W^ are 
shown at[T7]for the same set of lattice spacings. Inspect- 
ing these two figures, we do not see the strong volume 
dependence of the mp from the magnetic field. 



In addition and for the clarity we present at [18] the 
extrapolated value of p meson mass versus the squared 
magnetic field. The points are connected by splines. The 
lattice volume artefacts are not strong, but we are re- 
stricted by not sufficiently small lattice spacing and can 
not explore the masses at very big magnetic fields. As 
we know the radius of the lowest Landau level 



1 



rL 



(21) 



We could suppose that lattice artefacts may appear at 
r^ — a, so the simple estimations give for the spacing 
a = 0.0998 fm the maximum value of the magnetic field 
eB = 3.9 GeF^ 2.9 GeV^ for the a = 0.1155 fm and 
2.0 GeV'^ for the 0.1338 fm lattice. 

But from Figure [THJ we see that lattice spacing arte- 
facts appear even earlier, at ~ 1.5 GeV for the 0.1338 fm 
lattice and ^ 2~ 2.5 GeV for the other lattices. Notice 
that for the same fields and lattice volumes the masses 
for the finner lattices spacing lie a litle bit higher on the 
plot than the masses for the coarser lattices. So we can 
say that the neutral p meson mass doesn't become zero 
for the fields eB < 2 GeV^ and smoothly decreases. The 
larger fields demands bigger and finner lattices. 
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FIG. 16: The nip^^ in depence of bare lattice quark mass 
for different lattice spacings, magnetic fields and 14'' lattice 
volume. The extrapolated values of masses are depicted by 
stars and correspond to the physical quark mass m^'' = 3 
MeV. 



VIII. CONCLUSIONS 

In this work we explore the masses of tt, p and A 
mesons in the background of strong magnetic field of the 
hadronic scale. We observe that the masses extracted 
from different Lorentz components distinguish. The mass 
components extracted from the correlators of currents in 
the third direction along the field decrease with the value 
of the magnetic field, while the masses calculated from 
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FIG. 17: The same that at Figure [H but for 16 lattice 
volume. 
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FIG. 18: The nip.,., versus magnetic field for 14* lattice volume 
and a = 0.1155 fm lattice spacing. The values of masses are 
extrapolated to the physical quark mass rriq'^ = 3 MeV. 



the currents in the directions perpendicular to the exter- 
nal magnetic field slowly increase for p meson and stay 
almost constant for the axial meson. Indeed the symme- 
try in the perpendicular directions leads to the an ~ 022 



and pii ~ P22- We consider this phenomena to be the 
result of the anisotropy, which the strong magnetic field 
creates. We notice that our correlators give the masses 
of the neutral (not charged) p mesons and the neutral 
axial mesons. The phenomena of superconductivity cor- 
responds to the condensation of charged p - mesons [2^ . 
But the superconductivity of the charged mesons, which 
takes place at high values of field B, also have to induce 
the phenomena of superfluidity, this effects was consid- 
ered in QCD. The decreasing of p meson mass corre- 
sponds to the approaching of the superfluidity, but it 
seems the choice of the correlators doesn't allow to see 
clean so subtle effects, because of the following reasons. 
At first the averaging over the time slices strongly sup- 
presses the superconductivity and superfluidity signals 
in our correlators, because of an inhomogeneity of the 
ground state perpendicular to the magnetic field. With 
the growth of the field value these states become more 
wavy and the corresponding averaging should kill such 
states faster. However for the not very big fields such 
suppression by averaging can be irrelevant and we see 
the decreasing of the mass along the magnetic field. So, 
we conclude that the observed behaviour can be an evi- 
dence of the presence of superconducting and superfluid 
phases at the high values of the magnetic field B [30| • 
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